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ABSTRACT 


Stationary hot and cool particle distributions in the auroral magneto- 
sphere are modelled using adiabatic assumptions of particle notion in the 
presence of broad-scale electrostatic potential structure. The study has 
identified geometrical restrictions on the type of broadscale potential struc- 
ture which can be supported by a multispecies plasma having specified sources 
and energies. Without energization of cool thermal ionospheric electrons, a 
substantial parallel potential drop cannot be supported down to altitudes of 
2000 km or less. Observed upward directed field-aligned currents mist be 
closed by return currents along field lines which support little net potential 
drop. In such regions the plasma density appears significantly enhanced. 
Model details agree well with recent broad-scale implications of satellite 
observations. 
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1 . INTRODUCTION 


Over the past decade a quantitative description of the auroral quiet-time 
magnetosphere has been compiled from a large data base of satellite and 
ground-based observations (see Chiu et al., 1983; Burch et al., 1983; Mizera 
et al., 1981; Robinson et al. , 1982; and many others). Large scale features 
are represented schematically in Fig. 1. Detailed analysis of data from such 
sources has resulted in significant progress toward the identification and 
interpretation of plasma boundaries, particle energy and flux characteristics, 
electrostatic wave generation, and indeed, the entire spectrum of magneto- 
spheric phenomena. 

The carriers of the large-scale Blrkeland current linking the ionosphere 
and the magnetosphere have typically been identified as precipitating auroral 
electrons in the upward-current region and escaping thermal electrons in the 
downward-current regions (Anderson and Vondrak, 1975; Klumpar, 1979). Both 
species coexist, however, and often low-energy (< 50 eV) electrons can balance 
about half the primary upward current (Maier et. al., 1980). Narrowly colli- 
mated electron beams are often seen at the edges of auroral arcs (Anderson, 
1979), at energies below the parallel potential drop within lnverted-V's 
(Burch et al., 1979, Lin and Hoffman, 1979); and in the cusp (Zanetti et. al, 
1981). Upward-flowing and counter-streaming (100 eV) electron beams have 
recently been described by Lin et. al., (1982, 1984) and by Sharp et. al. 
(1980). There is some evidence to suggest that imbalances in the 
counterstreaming fluxes provide the effective field-aligned current flow in 
these cases, and that the width in pitch angle of the upward electron flux 
increases with increasing electron energy. Burch et. al. (1983) show that the 
pitch angle spread is related to the altitude of the top of a parallel 
potential drop below the spacecraft. 
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Figure 1 . A Schematic Diagram Summarizing Typical Observations and 

Interpretations Along Adjacent Evening Flux Tubes of Auroral 
Field Lines for Upward and Downward (Return) Current, 
Respectively. 


2 



Johnston and Winnitigham (1982) and Klumpar and Heikkila (1982) also 
( reported upward electron beams, the latter observations imbedded in a very 

broad region of downward Birkeland current. Collin et. al. (1982) showed that 
the occurrence frequency of upward beams is roughly eight times higher above 
6000 km than below 3000 km; suggesting a sporadic (nonadiabatic) low-altitude 
acceleration mechanism while a steady mechanism (perhaps a broad-scale poten- 
tial drop) dominates near 6000 km. 

Similarly, observations of ions in upward-current regions have identified 
conic (transversely-accelerated) ion distributions with energies of greater 
than 500 eV at a few thousand km altitude (Sharp et. al. 1977) , while upward 
ion beams are seen primarily above 6000 km (Gorney et. al. , 1981). Lower 
energy conic ion distributions were observed to originate at lower altitudes 
irrespective of local time or parallel current direction (Gorney et. al. , 
^ 1981; Whalen et. al. , 1978; Klumpar, 1979; Horwitz et. al. , 1981, 1982; Collin 

et. al. 1982). It is believed that the conic distributions are indicative of 
localized (1000 km thick) regions of non-adiabatic ion heating, while the 
upward ion beams reflect electrostatic acceleration by upward-directed paral- 
lel electric fields. Both mechanisms may act simultaneously (Klumpar et al. , 
1984). 

Due to the highly variable nature of the medium, individual data sets are 
often analyzed in terms of the local environment only. Nevertheless, it has 
been well established (Leraaire and Scherer, 1973, 1983; Mizera et al. , 1981; 
Chiu et al., 1983) that, for undisturbed conditions, the shapes of observed 
particle distributions are readily recovered by simple kinetic models which 
conserve the first two adiabatic invariants along auroral field lines. 
Kinetic theory qualitatively reproduces observed loss-cone, beam and conic 
features of individual particle species above 1000 km altitude. 
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In this paper, observed features of particle distribution functions along 
auroral magnetospheric field lines are reproduced by a two-dimensional sta- 
tionary model. The model is initialized by a representation of cross-field 
potential structure suggested by Chiu et al. (1981) as a solution allowing 
current closure within a finite domain of dosed field lines. Simple bi- 
Maxwellian distribution functions (as in the one-dimensional study of Chiu and 
Schulz, 1978) have been adiabatically mapped along field lines, taking into 
account acceleration due to monotonic parallel potential drops. The distribu- 
tion functions are then integrated in velocity space to provide particle 
densities as functions of magnetic induction and electrostatic potential for 
eight species of ionized particles. Shooting techniques are employed to 
identify smooth solutions for the electrostatic potential which satisfy quasi- 
neutrality everywhere. 

Section 2 describes the model and techniques for providing satisfactory 
solutions for the auroral potential structure. Density variations for eight 
particle populations are limited by adiabatic invariant boundaries in velocity 
space, and by effective particle temperatures, while they appear to be largely 
insensitive to our choice of bi-Maxwellian distributions. While such a simple 
direct model cannot accurately represent details of particle heating processes 
and pitch angle diffusion, these effects are induced through increased parti- 
cle temperatures and the existence of trapped particles. Hence, recognizing 
that significant electric fields are responsible for localized nonadiabatic 
heating, we employ model representations of ionospheric species which are 
assumed to behave adiabatically after experiencing heating in a low altitude 
acceleration region. 

Section 3 focuses on new interpretations based on these solutions. Our 
two-dimensional model results predict density cavities along field lines which 
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I support large potential drops (as seen by Calvert, 1981), and conversely, 

density enhancements along field lines which have no potential drop (i.e., at 
the edges of inverted V's as seen by Anderson, 1979). In regions of upward 
current, quasineutrality demands that any parallel potential drop along the 
lowest few thousand kilometers of the field line must be of the same order of 
magnitude as the temperature of the cool electron species. We show, in fact, 
that both the breadth and shape of the low altitude electrostatic potential 
structure appear to depend strongly on the energy of a population of cool 
electrons of ionospheric origin. In general, we expect that the potential 
drop along the lower part of a field line is never much greater than the 
temperature of the ionospheric species retarded by that potential. 

The concluding section provides interpretations of the distinct features 

of the model in terms of satellite observations, and discusses implications of 
v 

the model in regions of sharp density gradients at the edges of observed 
Inver ted-V events. The model will be shown to represent well the salient 
features of the data, Including prominent features predicted by the model and 
only recently experimentally confirmed. These features Include the identifi- 
cation of cool electrons of ionospheric origin as the current carriers in the 
low altitude return current region, enhancement in the density of most species 
along field lines which do not support a significant potential drop, and a 
strong correlation between the effective temperatures of particle species 
originating in the ionosphere and the depth of penetration of V-shaped poten- 
tial structures into the ionosphere. 
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2. TWO-DIMENSIONAL ELECTROSTATIC MODEL 


Aa in earlier studies (e.g., Chiu and Schulz (1978)], the particle motion 
is constrained in the magnetosphere by the conservation of energy ana magnetic 
moment along field lines: 


W = | (v , 2 + Vj 2 ) + q|e| * (1) 

M = m Vj 2 /^ B) , ( 2 ) 

where W is the total energy, m is the particle mass, q(" + or - 1) is the 
particle charge in units of the electron charge |e|. are particle 

velocities parallel and perpendicular to the magnetic field, respectively 
(Vg > 0 is downward). 41 and B are the electrostatic potential and magnetic 
field magnitude, respectively. 

The simplest particle distribution functions are characterized at a 
convenient point of origin as bi-Maxwellians: 

f a exp (W| /Tj + VV 


or 

2 v 2 

f a exp [j (—— + + 9 |e| £-) • (3) 

II 1 n 

Once the effective temperatures and electrostatic potential are determined at 
one point, which we designate as the '’source" of the particle distribution, 
then the adiabatic invariants associated with the conservation of energy and 
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magnetic moment c<« . ' ately determine the particle distribution at every other 
point along the magnetic field line — as long as the electrostatic potential 
is known everywhere. 

Boundary conditions ore of great importance. Assuming conjugate-point 
symmetry along field lines, we define the upper boundary of the system to be 
the magnetospheric equator, the source region for hot plasma sheet protons and 
electrons. Particles traveling downward which reach the lower boundary are 
presumed to be lost from the system. The lower boundary at the base of the 
field lines is assumed to be an ionospheric sink (or source) located high 
enough that particle collisions may be neglected. We have chosen 1000 km 
altitude as the altitude in our model at which e*».bl polar diffusion is no 
longer dominant and magnetospheric quasineu?rali fy? limits particle flux. 

The region below 1000 km is complex due to colllslonal effects. To avoid 
complicating our representation of the magnetospheric plasma, we chose this 
altitude as a lower boundary characterized by an ionospheric source which 
models the low altitude plasma. Cool ions below 1000 km altitude, which are 
gravitationally and colllslonally bound, slow the escape of electrons up the 
field lines. An effective upward-directed electric field component is 
generated by ambipolar diffusion in this region. Note that this region 
effectively restricts most of the ionospheric particles to the ionosphere 
because only the most energetic can escape. This also answers a complaint 
voiced by early opponents to the idea of higher altitude parallel electric 
fields directed downward; scientists who felt that such fields would evacuate 
the ionosphere. The ionosphere is preserved because downward electric fields 
do not extend through the ionosphere to the earth's surface. (The reversed 
potential drops will be examined in detail in work to be presented 
elsewhere.) Gradual heating of the ion and electron populations occurs as 
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electrons are energized by the precipitating flux, and ions are energized by 
perpendicular field turbulence in the form of wave/particle interactions. 

In this paper we sidestep the question of low-altitude local heating 
mechanisms to address the question of the accumulated affect of this heating 
on the high altitude distributions. How hot oust the ionospheric species be 
in order to provide a quasineutral plasma at high altitude with the 
distribution function boundaries seen by satellites in these regions? Our 
cooler ionospheric distributions in this report represent only those particles 
which could escape the ionosphere, that 1 b, superthermal particles. 

The bi-Max>wllian distribution functions are further restricted by an 
imposed accessibility criterion which stares that the parallel velocity as 
defined by (1) must be real at the source, and at all positions along a field 
line between the source and any position, s, where the density is nonzero. 
This states that a particle cannot pass through a region where the probability 
of particle occurrence is zero. Note that even abrupt monotonic changes in 
the electric field are acceptable as long as the particle is accessible to the 
point of Interest, 

Particle v&? ocAty-space distribution functions are restricted to specific 
regions of velocity-space by imposed boundary conditions and adiabatic invari- 
ants. The adiabatic Invariants relate the velocity distribution at a boundary 
to that at any position s ( see Chiu and Schulz, 1978; recall that ions 
require an extra term to include the effect of gravity): 

v u 2 “ v iis 2 + U “ (V b s)K 8 2 + U s ^ £ ) (5) 
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Hence, the boundary limitations (s * 1,0 ) may be plotted in velocity space for 

any particle at position s, where l refers to the lower boundary and 0 to the 

magnetospheric equator* This is illuserated in Fig, 2 for hot magnetospheric 

electrons accelerated by a 2 keV drop in the potential along the length of a 

field line. The hot electrons are restricted to the region outside the 
2 

ellipse VgQ “0; that is, all must have access to their point of origin at s 

2 

" 0. The hyperbola, v^ "0, separates sections of velocity space occupied 
by particles of magnetospheric origin which mirror, from that occupied by 

particles which precipitate. In applying these boundary conditions, particles 
which mirror before they reach the lower boundary (v. < 0) are 

R X 

distinguished from those which precipitate through the lower boundary and are 
lost (v 58 > 0, v u 2 > 0). 

As shewn in Fig. 2a, nsost hot elections (of several keV) near the origin 
are mirroring electrons (cross-hatched area). When a parallel potential drop 
of 2 keV is assumed along the field line, nuch more of velocity space is 

accessible near the top of the field line (Fig, 2a) than at lower altitudes 
(Fig. 2b). A large fraction of those electrons which are able to approach the 
bottom boundary (Fig. 2c) also have a high downward (v^> 0) parallel component 
of velocity and will precipitate. Electrons which pass down the field lines 
through a sizeable pot. '.tial drop (upward E ) are all accelerated to higher 
speeds, leaving the central ellipse of the distribution function empty (see 
Fig. 2c); describing partial ring distributions. The ellipse is not filled by 
electrons streaming out of the ionosphere because upward motion of cooler 

electrons is effectively curtailed by the potential drop. The realization of 
this fact led researchers to interpret observed inverted-V events as distinct 
signatures of parallel potential drops. If, however, mechanisms exist to 
energize these cool electrons to 50-2000 eV, electron conics are seen 

travelling obliquely up the field lines [Menietti and Burch, 1985], 
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Figure 2. Regions of Velocity Space Occupied by Electrons of 

Magneto spheric Origin, Which Either Mirror Back Up Diverging 
Field Lines or Precipitate Toward the Lower Ionosphere. A 
potential drop of 2 keV is assumed along the field line. 
Regions are shown for particle distributions a) near the 
source of particles (magnetospheric equator) ; b) at an 
Intermediate position; and c) near the base of the field 
line (assumed to be at altitude of about 1000 km) . 
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In observed particle velocity-space distributions non-adiabatic processes 
(e.g. t pitch-angle diffusion) smooth the edges of the boundaries in velocity 
space. Nevertheless, it is possible to distinguish contributions to different 
regions, bounded by these curves as shown in Fig. 3a for ions, and Fig. 3b for 
electrons. If downward-directed electric fields exist, the boundary curves 
for ions and electrons are reversed as shown in Fig. 4. Typical distributions 
sample mixed populations of both hot and cool particles; the hot particles 
originating near the plasma sheet while the cooler particles have their source 
in the ionosphere. Assuming that the particles conserve adiabatic invariants, 
estimates of the potential drop necessary to produce observed distribution 
functions are routinely made from satellite data (Mizera et. al. , 1981; Chiu 
et, al. , 1982). 

Here we extend the theory to describe particle distributions along the 
entire length of magnetic field lines across the auroral zone. The densities 
of all representative particle populations are obtained by integrating over 
adiabatically allowed regions of velocity space and applying shooting tech- 
niques to solve for values of the electrostatic potential which satisfy 
quasineutrality among the particle densities at every point in space. These 
models contain both regions of upward current and regions of return current in 
order to preserve charge balance between the magnetosphere and ionosphere. 

Our technique for achieving self-consistent solutions requires initial 
estimates of the two-dimensional potential structure. The structure along the 
top boundary (magnetospheric equator) is assumed to be a smoothly-varying 
function of cross-field distance. By choosing a simple potential well, we 
have been able to confirm analytically our numerical solutions for the 
potential at the lower boundary which satisfy current continuity within the 
system (Chiu et. al. , 1981). Conceivably, any potential structure could be 
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Figure 3. Boundaries in Velocity Space for Various Particle 

Populations at an Intermediate Position Along a Field Line 
in a Region of Upward Current. The species are either a) 
ions, or b) electrons designated by M if they are assumed 
to originate in the plasmasheet (magnetosphere) , and by I 
if the assumed origin is the ionosphere. B designates 
backscattered electrons, and T designates electrons 
trapped between coverging magnetic field lines and a 
parallel potential barrier. The dashed diagonal lines 
are asymptotes Vj = “ l] 1 ' v^ t ^ e B °H-d 

hyperbolic boundary ■ 0. 





e (<t> ~ 4 >q) \ " 

rrr V®~/_ 


1/2 


[From Chiu and Schulz, 1978]. The scale size for ion 
velocities is ten times the electron velocity scale. 
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Figure 4. 


Boundaries in Velocity Space as per Figure 3, However, 
for a Reversed Parallel Potential Drop, Hence a Return 
Current Region. The region designated T indicates 
that trapped ions may exist at some positions along a 
field line. The scale size for ion velocities is one 
hundred times the electron velocity scale of Figure 3. 


v^ and Vj^ have the form given in Figure 3 with 
m_ replaced by m + . 
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applied at the upper boundary; our simple choice is based on an understanding 
of typical length scales across auroral field lines near the equator. For the 
ionospheric potential <J» , the height-integrated current continuity equations 

f* 

reduces to 


h • (E P ’i ♦*> " - J . 

where is defined as negative for downgoing electrons and is the height- 
integrated Pederson conductivity. Kinetic models and observations imply 
that Jg is proportional to the field-aligned potential difference between the 
magnetospherlc equator and the ionosphere. By estimating the proportionality 
constant, one can express the current conservation equation as a Poisson-like 
equation which allows us to estimate the potential structure at the lower 
boundary (as in Chiu et al. , 1981). An interpolation linear in B then 

provides an initial estimate of the electrostatic potential everywhere between 
the boundaries. This interpolation procedure is chosen because the kinetic- 
theory proportionality of potential to J along a field line, yields a 
potential proportional to B, which leads to a current conserved without 
transverse components. Recent experimental observations by both DE spacecraft 
along the same field line confirm this proportionality (Weimar, 1985]. 

Once the potential is estimated everywhere, the particle density func- 
tions are obtained by integrating the distribution functions over allowed 
regions of velocity space for all representative particle populations. These 
density functions depend only on magnetic field, electrostatic potential, 
effective particle temperatures (parallel and perpendicular to the magnetic 
field), and on a constant multiplicative coefficient defined by boundary 

conditions. At the upper boundary a balance is assumed between the densities 
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of oppositely-charged particles of plasma sheet origin. At the lower bound- 
ary, a similar condition forces charge neutrality among particles of iono- 
spheric origin. Initial estimates of the densities are made for a variety of 
particle temperatures within the energy ranges suggested by satellite data. 
The total charge density is calculated on a coarse grid between the bound- 
aries, and parameter sets describing the particle temperatures and density 
coefficients are chosen for which quasineutral solutions appear likely at all 
spatial positions. 

Quasineutral solutions are quite adequate for the range of parameters of 
interest in the evening magnetosphere, as long as variations on spatial scales 
comparable to the ion Larmor radius are not of concern [Chiu et al. (1983) J. 
Sharp variations in the electric field are not of interest because they would 
invalidate the assumptions of adiabaticity inherent in our model. The 
smoothly- varying large-scale potential structure described by our model is 
representative of the net charge separation and E associated with inverted-V 
structures, where the actual net charge density on auroral field lines is very 
small. If a quasineutral solution could not be obtained everywhere across 
such a large-scale system, then it would also be unlikely that the derived 
particle densities could adequately describe a steady, quiet magnetosphere. 
Reasonable solutions are In fact Identified, particle densities are calculated 
at each point on a fine grid, and shooting techniques are used to refine the 
value of the potential until quasineutrality is achieved at each grid point to 
within an absolute error, 2|n + - n_|/(n + + n_), of 0.001 or less, 

A few statements must be included here to address the uniqueness of the 
solutions and the ease in obtaining them. We have found that once species* 
temperatures and coefficients are defined and boundary conditions are applied, 

then there will be either one physically meaningful solution, or none. Our 
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system of constraints Is too rigid to allow an unreasonable solution. It is 
feasible, for example, that shooting techniques could converge on a value of 
the potential inappropriate for low altitudes, in an attempt to balance, say, 
cool ions with trapped electrons rather than with cool electrons when the 
latter are too cold to be accessible to these heights. Invariably, in such 
cases the technique fails miserably near the lower boundary, and further 
efforts must be extended to find appropriate cool electron temperatures. A 
single solution at any point is not necessarily unique; however, a smooth 
consistent solution everywhere may often be obtained by predisposing the 
shooting technique to search for a solution near the value computed at the 
same altitude for an adjacent field line. Any discontinuous solutions are 
rejected as Inconsistent with adiabatic invariant assumptions. The parameters 
listed in Table 1 are representative of typical solutions which satisfy 
quasineutrality everywhere. The evening magnetosphere may adopt one of a 
number of steady configurations, differing in representative temperatures of 
different species, and in the shape of net potential structure. 
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Table 1. Parameters Defining the Particle Species for Six 
QuaBineutral Solutions 


Particle 

Example j 

Species/ 

Parameter 

1 

2 

3 

4 

5 

6 

"I- T ll 

1.2 

40 

1 

2 

.4 

.8 

T 1 

1.2 

40 

2 

2 

.4 

.8 

C(xl0“ 3 *) 

76 

13 

41 

59 

45 

45 

n H+ T | 

0.5 

0.6 

0.7 

5 

.03 

.03 

T 1 

10 

20 

20 

5 

.03 

.03 

C(xlO“ 35 ) 

3.9 

2.1 

2.3 

25 

3,4xl0 8 

3.4xl0 8 [ 


1 

3 

5 

10 

.15 

.35 

T 1 

15 

20 

50 

10 

.15 

.35 

C(xC^) 

30 

71 

34 

23 

4.2 

4.2 

n B- T 

16 

400 

400 

400 

300 

300 

C 

0.34 

0.34 

0.34 

0.34 

0.34 

0.34 

n>p_ T 

1500 

1000 

1000 

1500 

400 

400 

C 

61 

40 

32 

83 

64 

64 

Upward Current 
Solution 

X 

X 

X 

X 



Return Current 
Solution 





X 

X 


Temperatures are in eV. Multiplicative coefficients « C a in cc 
depend on distribution function form. n M+ is isotropic with T-6 
keV. n M _ is anisotropic in the upward current region with T «* 
1.5 keV and T^ - 3 keV at the central field line; but isotropic^in 
the return current region with T-2.5 keV, (see text). C,. = 1.1 x 
10“ 33 ; C M _ » 2.8 x 10“ 38 . 


18 




3. SMOOTH SOLUTIONS - SPECIES DENSITIES AND ELECTROSTATIC POTENTIAL 

We will consider in detail several sets of quasineutral solutions. The 

particle parameters differ for regions of upward or return current but are 

otherwise constant. Table 1 lists parameters describing quasineutral 

multicomponent plasma species under various magnetospherlc conditions. The 

effective temperatures quoted describe a bi-Maxvellian distribution for each 

species only at its source. Further along the field line adiabatic invariants 

may Impose a highly assymetric distribution. Two sets of parameters are used 

to define a standard magnetosphere, the choice o i upward and return current 

parameters determined by the requirement that the total upward current balance 

the total downward current of the complete system. (The current density is 
2 

2.6 pA/m at the base of the field line carrying the strongest upward 
current. ) 

In Fig. 5, examples 1 and 5 from Table 1 for the two current regions are 
compared with examples 2 and 6. Both electrostatic potential and total plasma 
density are computed on a representative 12 x 46 point grid in a meridional 
plane across the auroral zone. The model calculations employ dipolar geometry, 
however, for clearer illustration we have plotted the results employing the 

familiar transformation x * Xg (Bq/2B^)^^, compressing divergent field lines 
to parallel vertical lines with the magnetospherlc equator at the top boundary 
and 1000 km ionospheric altitude at the bottom. The cross-field distance 

scale refers to distances between field lines at the lower boundary (the 

baropause at 1000 km altitude). The transformation acts as an image 

enhancement, readily identifying gradients in the potential or density with 
respect to the magnetic field geometry. The shaded central portion of each 
plot identifies a region of upward current surrounded by unshaded symmetric 
return current regions. Note that the physics does not specify whether return 
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Figure 5. Computational Solutions for a) the Electrostatic Potential 
(kev) and b) Total Electron Density (/cc) Along a 
Geomagnetic Meridional Plane. Magnetic field lines are 
plotted as vertical lines in this projection. As 
labelled, the upward current regions (shaded) arc 
described by Table 1 examples 1 and 2. Similarly, the 
unshaded regions are symmetric return current regions 
described in Table 1 by examples 5 and 6. 
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currents should be north or Bouth of an upward current region. Certainly both 
are seen, so we include both in our representations. 

Fig. 5 illustrates the kind of broad-scale variations to be expected in 
satellite observations. Narrow potential structures which become increasingly 
narrow at lower altitudes may be seen in conjunction with localized plasma 
density cavities which become smaller at lower altitudes. Broader and deeper 
potential structure correlates with more extensive plasma depletions, bounded 
by very narrow density enhancements. The character of specific solutions 
stems directly from the response of individual particle species to the global 
electric and magnetic field structure. 

Fig. 6 illustrates the density variation of three separate electron 
populations. The hot electrons of magnetospheric origin are isotropic (2.5 
keV) in the return current region, but are allowed to become gradually 
anisotropic toward the central field line where (T ([ , T^) ■ (1.5, 3) keV.* This 
feature has been incorporated to match interpretations of observed 

velocity-space anisotropies, but it has only a minor effect on the density 
structure of the hot electrons, and practically no effect on quasineutrality 
considerations. These hot electrons are most dense at their origin, the top 
boundary of Fig. 6a. The electrons are least dense in the central region 
where they are accelerated downward by a potential drop of 2600 eV, and many 
are lost due to precipitation through the lower boundary. Reduced 
precipitation along equipotential field lines near the outer edges of the 
upward current region produce a density maximum there. The return current 
region is assumed to support a small parallel potential drop which is assumed 
by us to be steady over hundreds of kilometers, producing a constant density 
versus latitudinal distance. The return current potential drop is 

characteristic of an extended ambip'j potential, and is sufficient to 
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Figure 6. Solutions for the Density (/cc) of electrons a) of 

Magnet ospherlc Origin, b) Trapped Along Flux Tubes, c) of 
Ionospheric Origin. Numbers In brackets denote Maxwellian 
temperatures (T , T ) at the species origin In the upward 
or return current regions as Indicated. The format Is as 
given In Figure 5. 
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balance the upward current when the sum of the small currents carried by each 
field line is evaluated* 

Trapped electrons (Fig. 6b) occupying the region of velocity space 
labelled T in Fig. 3, arc restricted to the upward current region and are most 
dense near the lower boundary. Double maxima occur on either side of the 
central field line at several thousand kilometers altitude. The central 
minimum between the maxima in the density contours of Fig. 6b occurs because 
the central field lines carry a more significant potential drop to very low 
altitudes, hence further restricting the allowed region of velocity space for 
trapped particles along these i... 3 ld lines. Technically, there is no mechanism 
to produce trapped particles in a completely adiabatic system, but we know 
that diffusive processes tend to fill restricted areas of velocity space. We 
assume that the population exists, and that it is maintained by the gradual 
diffusion of electrons across boundaries in velocity space; a minor 
perturbation on an otherwise stationary adiabatic system. The model solutions 
for the trapped electron distribution functions suggest that these are 
basically electrons of ionospheric origin which have received significant 
transverse acceleration at low altitude and have diffused upward into a region 
where they are trapped adiabatically between a potential barrier at high 
altitude and the dipolar magnetic mirror at low altitude. 

Untrapped electrons of ionospheric origin are represented in Fig. 6c. 
These electrons are so dense at low altitudes that they essentially determine 
the net character of density gradients in the plasma as a whole. The numbers 
in brackets above the lower boundary are the effective bi-Maxwellian parallel 
and perpendicular temperatures (isotropic for these examples) defining this 
cool electron population at the lower boundary. In both examples given, the 
assumed temperatures are super thermal. The order-of-magnitude difference in 
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their temperatures accounts for the difference in the gradients of both 

plasma density and potential of Fig. 5. The maximum in this electron 
distribution is always at the outer edge of the upward current region along 
field lines which support no potential drop. 

In the upward current region where cool Iona are accelerated upward by 
the central parallel potential drop, it is necessary to have a mechanism to 
allow a sufficient number of cool electrons to follow the ions so that 
quasineutrality is satisfied between these dense, cool ionospheric populations 
at lew altitudes. Very cool electrons cannot cross a significant potential 
barrier, so the low altitude potential drop must be very small in the first 
example of Fig. 6c. Warmer electrons can penetrate further into a significant 
potential structure; or, to look at the phenomena another way, the potential 
structure can penetrate lower in the magnetosphere without retarding the 

ionospheric electrons too aeverely. In both examples, a density minimum in 
this population occurs along the field lines supporting the largest potential 
drop. Such density cavities have been observed by Calvert (1981) and others. 
They suggest that density cavities are correlated with observations of auroral 
kilometric radiation (AKR), possibly identifying source regions for electron 
cyclotron energization. In the examples of Fig. 5 the central density 
cavities do meet the threshold condition for AKR; that is, the plasma 

frequency is less than one fifth of the electron cyclotron frequency. Warm 
electron distributions within the density depletion appear conic-like as 

observed recently by Menietti and Burch (1985). If the potential drop is 

small, dipolar magnetic field geometry causes cooler electron distributions to 
appear to be upward-travelling beams or counter-streaming beams above 3000 - 
4000 km, agreeing with observations of Lin et al. (1984) and Sharp et al. 
(1980). 
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One might expect an important role to be played by another electron 
population energized at low altitudes through backscatterlng by hot 
precipitating electrons* Backscattered electrons are Included In the model, 
based on predictions of Prasad et al. (1983). These predictions suggest a 
maximum density of less than 0.5 particles per cc, that is, too small to 
significantly affect quasineutrallty. The backscattered electron density 
contours are not plotted In Fig. 6 because they are a very minor component. 
The density contours have similar shape to scaled trapped electron predictions 
without the central minimum. Cool ionospheric, trapped, and backscattered 
electrons all serve to neutralize the accelerated ionospheric ions streaming 
upward in the region of upward current. Distinguishing between these electron 
populations is difficult experimentally. The ionospheric electrons and 
backscattered electrons occupy the same region of phase space and are hence 
indistinguishable from each other on the scales of importance here. Our 
models Indicate only that warmer-than-thermal electrons are necessary to 
achieve a quasineutral state in the 1000 - 4000 km region. Our representation 
of the balance between backscattered and cooler electrons may be considered 
tentative, in liu of more quantitative data which is able to distinguish 
between them. It is clear, however, that if the potential drop extends to 
great depth in the ionosphere some nonadiabatlc process must exist which 
transfers sufficient energy to the cool electrons to allow them to reach 
adequate heights following the accelerated ions. 

The most important ion populations corresponding to the electrons of Fig. 
6 are illustrated in Fig. 7. The magnetospheric ions in the model are assumed 
to be isotropic at temperatures of 6 keV. These ions are only slightly 
retarded by the potential drop, the lowest density seen in the center of the 
upward current region. Trapped ions may also exist but only in the return 


25 


I) density of protons of MAGNETOSPHERE ORIGIN 




t» density of protons of ionosphefiic origin 



B 

(G) 



DISTANCE (x 100 km) 


C) DENSITY OF OXYGEN IONS OF IONOSPHERIC ORIGIN 


DISTANCE (x 100 km) 


Figure 7. Solutions for the Density (/cc) of Ions: a) Protons of 

Magnetcspherlc Origin, b) Protons of Ionospheric Origin, 
c) Oxygen Ions of Ionospheric Origin. The format is as 
given In Fig. 6. 
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current regions which support a reversed parallel potential drop* Even within 
a reverse potential drop it is unlikely that they have a major role in charge 
balance. Observations of low energy ions Indicate beam and conic shaped 
distributions rather than the shape predicted for trapped distributions in 
Fig. 4a. Although some argument can bo made that the observed lonB are the 
result of non-adlabatlc heating events at magnetospheric altitudes of 
1000-5000 km altitude, it 1 b clear that the bulk of the distributions seen do 
not correlate with trapped ion distributions along field lines above a 
localized heating region. The beam and conic ions appear to map essentially 
adiabatlcally as warmed distributions of ions that are not significantly 
retarded by any potential barrier! One suggestion that comes tr. mind is that 
the sluggish ions absorb energy from local wave fields (through resonant, 
non-adiabatic processes which tend to deplete the wave field). The warmed ions 
are subsequently channelled by the dipolar magnetic field geometry into conic 
distributions at low altitudes and beam distributions at higher altitudes. 
Note that movements of localized field structure could cause variations in the 
height of the heating region for ions, and indeed many such regions could 
feasibly exist. 

The resulting ion distributions can be modelled by recognizing that these 
effects tend to smooth the particle distributions. With minor distortions at 
very low altitude, the observed distributions are represented as somewhat warm 
populations of ions originating at low altitude. Both low-energy hydrogen and 
oxygen ions are assumed to exist. In the example of Fig. 7, highly 
anisotropic distributions were assumed in the central acceleration region with 
perpendicular temperatures exceeding parallel temperatures by factors of 7 to 
33 as shown. These ions appear as beams and conics, originating near the 
lower boundary (1000 km altitude), after presumably having undergone 
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non-adiabatic heating in some altitude range near this boundary. Recall, the 
non-adiabatic process is not represented here, just the result of the process, 
that is, the existence of warmed ions. 

The most distinct features of the ionospheric ion distributions shown in 
Fig. 7 are the sharp increases in the density along field lines which have no 
potential drop. Since these ions are not accelerated out of the region by an 
external electric field, they are more dense. Technically, there are more 
ions along these field lines because fewer are lost at the boundaries (more of 
velocity - space is assessible to them). In this model we consider the 
ionospheric source region at the bottom boundary to be a constant; hence, 
particles returning to the source do not supplement the Bource. The model 
thus predicts a steady-state distribution of ions with more ions attached to 
some field lines than others. (This same argument is seen to be valid also 
for the ionospheric electrons.) Other features of the model include a 
decrease in ion density along field lines supporting the largest potential 
drop. Both ions and electrons spend proportionately less time in such regions 
of strong electrostatic acceleration. Above about 2000 km the cool 
ionospheric electron population is excluded by strong fields in the central 
region, and the ionospheric ions and trapped electrons form a quasineutral 
balance. 

The altitude above which cool ionospheric electrons are negligible 
differs in the two examples shown because of details in the low altitude 
potential structure. In Fig. 8a latitudinal variations in the electrostatic 
potential are plotted at constant values of the magnetic field (ie. 
essentially constant altitude). The base of the field line has a magnetic 
field of 0.38 G, while the magnet os pheric equator is at 0.00053 G. Following 
a curve for a constant magnetic field of say 0.14 G (4000 km altitude) in 
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example 1/5, we see that the potential matches the lower boundary potential 

for all but the central half of the upward current region. That is, up to 

4000 km the outer half of the field lines in the upward current region are 
equipotentials — motion up the field lines is not Impeded for any cool 
species, and particle densities are higher than they would otherwise be. In 
example 2/6, there is a small fairly constant potential drop in the same 
region. The more energetic part of the 40 eV cool ionospheric electron 
population could crosB this small potential barrier. Quasineutrality is 
satisfied if there are not too many cool ions reaching this altitude. Fig. 6, 

7 and 8 show that both electron and ion densities are lower in this region for 

example 2/6 than for example 1/5. Small variations in the potential drop at 

low altitude have a strong effect on the distribution of cooler particles at 
higher altitudes. 

High altitude variations in the potential drop are more difficult to 
interpret. Fig. 8b illustrates variations In the potential drop from the top 
boundary as the auroral zone is crossed in the not tlvsoul.h direction. If a 
satellite maintained a constant altitude (approximately constant magnetic 
field) over several degrees latitude N-S distance, it would observe a fairly 
constant potential drop from the top boundary, throughout the central portion 
of the upward current region down to an altitude of about 4000 km (B a 
0.14G). This does not suggest that downgoing particles of magnetospherlc 
origin would contribute constant densities at a specific altitude in the 
central region. The structure near the top boundary in the graph of 
magnetospherlc electrons (recall Fig. 6) is caused by mirroring hot electrons 
which are also influenced by the variable potential drop belo„ this 
altitude. Nonetheless, one general rule which holds for all untrapped 
particle species is that lower density is found along field lines which 
support a significant potential drop. 
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4. IMPLICATIONS FOR PARTICLE HEATING 


Observations of warm ion populations (conics) directed obliquely upward 
(Gorney et. al. , 1981; and Klumpar, 1982) have helped to define the 

anisotropic ion distributions in our model. Recent observations have 
highlighted striking differences in the motion of different ion species 
[Hultqvis t, 1983; Horwitz, 1982). Differences in heating regions and heating 
rates are Important; however, we can readily show that a factor at least as 
important is the difference in adiabatic motion of different species within 
static fields. This is well illustrated by plotting the boundary restrictions 
for both oxygen (dotted line) and hydrogen (solid curve) in regions for which 
an upward/downward (Fig. 9a/b) parallel electric field component exists 
between the field line position (s) and the ionospheric boundary. In an 
upward current region, as in Fig. 9a, individual hydrogen ions (H) would be 
accelerated to greater velocities by the potential drop than would the oxygen 
ions (0). In addition, oxygen ions must overcome a gravitational barrier 
sixteen times as great as for hydrogen ions, so the net accelerating potential 
is greater for the hydrogen ions. As a result, the hydrogen ions spend less 
time in the vicinity of position s, and would he less dense there (all other 
parameters being equal). Although the ratio of the densities of oxygen to 
hydrogen ions at the lower boundary is taken to be a half, the ratio increases 
to a maximum of 1.1 higher up the central field line over the arc. If ion 
heating occurred at low altitude, it is clear that the oxygen ions would spend 
much more time in this region than the hydrogen ions simply due to their 
slower parallel velocities. Local heating would thus allow an even greater 
proportion of oxygen to escape further along the field lines. Note that in 
the examples we consider, the accelerating potential drop extends down to the 
lower boundary over the central field line. If conditions were very quiet, 
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al 0 V H b) 0 


Figure 9. Representations of 0 4 and H 4 Distributions near Their 

Ionospheric Origin in a) an Upward Current Region, b) a Return 
Current Region with Reversed Potential Drop. The ratio of 0 4 
to H 4 varies with altitude, parallel potential drop and the 
anisotropy (T /T ( ) of the density distributions at the 
origin. The ion velocity scale is consistent with Figure 3. 
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the potential drop would not extend as deep, and It is likely that a much less 
significant oxygen ion component would petie crate into the magnetosphere. 
Prognoz observations described by Hultqvist (1983) indicate that the oxygen to 
hydrogen ratio at geosynchronous altitudes may be a factor of ten or so during 
disturbed timeB. Horwitz (1982) also observes a stronger oxygen component at 
1400 km for very disturbed cases only. Most often the potential drop does not 
extend deeper than 5000 km (Gorney et al, , 1981) and the ratio is less than 
one. It is of interest that even in disturbed conditions the energy in the 
observed ojqrgen and hydrogen distributions (at geosynchronous altitudes) is 
about the same for both species; indicating that the dominant source of free 
energy is probably electrostatic in nature. 
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5. DISCUSSION 


In Chls paper, we have described models of electrostatic potential dis~ 
tributions consistent with current continuity and change balance, for particle 
populations governed by 'adiabatic invariants and quasineutrality in the 
magnetosphere. The models reproduce the general characteristics of quiet-time 
satellite observations over evening auroral arcs, both in regions of upward 
and return current. Estimates of the current density show that the upward 
current density above 1000 km is carried by a population of hot electrons of 
magnetospheric origin. (Although the trapped electrons may carry appreciable 
current in either direction, the net effect is zero.) Throughout the return 
current region, the current is carried by the cooler electrons of ionospheric 
origin. We also infer from the models that both the latitudinal extent of the 
electrostatic potential structure, aB well as the depth of penetration of this 
accelerating potential into the ionosphere are strongly influenced by the 
effective energy of the cool electron population. Our modelling results 
conclude that broad potential structure extending low along the field lines 
can only describe a stationary qujsineutral solution if the representative 
ionospheric electron population is assumed to be suprathermal (10 - 100 eV). 
Hence, signature of strong parallel electric field components seen to extend 
down to lower altitudes, c> oild be interpreted as a sign that non-adiabatic 
heating of the cooler electrons had taken place. Thick regions of enhanced 
density at the outer edges of inverted-V's may indicate a cooler ionospheric 
electron distribution. In this case quasineutrality demands that the field- 
aligned potential drop be suppressed at low altitudes (allowing some cool 
electrons up the field lines to balance ions of ionospheric origin). This 
also tends to reduce the crosB-field extent of the potential structure, 
allowing the magnetic field lines to appear as approximate equipotentials up 
to great altitudes. 
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